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According to general relativity, a spinning body of mass M and angular momentum S, like a star 
or a planet, generates a gravitomagnetic field which induces, among other phenomena, also the 
Lense-Thirring effect, i.e. secular precessions of the path of a test particle orbiting it. Direct and 
indisputable tests of such a relativistic prediction are still missing. We discuss some performed 
attempts to measure it in the gravitational fields of several bodies in the Solar System with natural 
and artificial objects. The focus is on the realistic evaluation of the impact of some competing 
classical forces regarded as sources of systematic uncertainties degrading the total accuracy ob- 
tainable. In the case of the test performed with the LAGEOS and LAGEOS II Earth's satellites 
one of the major sources of systematic uncertainty is the imperfect knowledge of the even zonal 
harmonic coefficients of the multipolar expansion of the Newtonian part of the terrestrial gravita- 
tional potential. The lingering uncertainty in some of them makes the total error in such a test as 
large as 15 — 30%, contrary to more optimistic evaluations (« 5 — 10%). Some different strategies 
to extract the relativistic signature from the data of the LAGEOS satellites which may comple- 
ment and robustly corroborate the so far implemented tests are suggested. We critically discuss 
the possibilities that the LARES satellite, to be launched at the end of 2009, will realistically al- 
low to measure the Lense-Thirring effect with a ss 1% accuracy. Since it will orbit at much lower 
altitude than LAGEOS and LAGEOS II, much more even zonals will have to be accounted for; 
according to the present-day uncertainty in them, calculations performed with standard geodetic 
techniques show that their impact may be orders of magnitude larger than the expected accuracy 
level. Recent progresses in the planetary orbit determination make the perspective of reliably 
measuring the Sun's Lense-Thirring effect a realistic possibility. Presently, the magnitude of the 
gravitomagnetic perihelion precessions of the inner planets is about of the same order of magni- 
tude of the present-day uncertainties in determining the secular perihelion precessions from the 
observations of the rocky planets. Moreover, the predicted Lense-Thirring effect for all of them 
are in agreement with the estimated corrections to the standard Newtonian/Einsteinian perihelion 
precessions. Finally, we discuss a recent interpretation of some orbital data of the Mars Global 
Surveyor spacecraft in terms of the gravitomagnetic field of Mars which recently raised a debate. 
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1. Introduction 

In the weak-field and slow motion approximation, the Einstein field equations of the gen- 
eral theory of relativity, which is a highly non-linear Lorentz-covariant theory of gravitation, get 
linearized resembling to the Maxwellian equations of electromagntism. As a consequence, a grav- 
itomagnetic field B g , induced by the off-diagonal components go;, i = 1,2,3 of the space-time met- 
ric tensor related to the mass-energy currents of the source of the gravitational field, arises [j|]. 
The gravitomagnetic field affects orbiting test particles, precessing gyroscopes, moving clocks and 
atoms and propagating electromagnetic waves |fj, BJ|. Perhaps, the most famous gravitomagnetic 
effects are the spin-spin precession of the axis of a gyroscope [|, |] and the spin-orbit Lense- 
Thirring 1 precessions [^] of the orbit of a test particle, both occurring in the field of a central slowly 
rotating mass like, e.g., our planet. For an analysis of the ambiguity arising from the use of the ter- 
minology frame-dragging, de Sitter precession 2 and Lense-Thirring precession, in contrast to the 



unambiguous reference to spin-orbit and spin-spin precessions, see Ref. Q10|]. Direct, undisputable 
measurements of such fundamental predictions of general relativity are not yet available. 

The measurement of the gyroscope precession in the Earth's gravitational field has been the 
goal of the dedicated space-based GP-B mission jTT[ |l2| ] launched in 2004 and carrying onboard 
four superconducting gyroscopes; its data analysis is still ongoing. The target accuracy was orig- 
inally 1%, but it is still unclear if the GP-B team will succeed in reaching such a goal because of 
some unmodelled effects affecting the gyroscopes: 1) a time variation in the polhode motion of the 
gyroscopes and 2) very large classical misalignment torques on the gyroscopes. For early results, 



see Ref. [13]; see also on the WEB http://einstein.stanford.edu/. 

In this Lecture we will focus on the attempts to measure of the Lense-Thirring effect in the 
gravitational fields of the Earth, the Sun and Mars. Far from a localized rotating body with angular 
momentum 5 the gravitomagnetic field can be written as 



Bg cr 3 



5-3 {S-f )r 



(1.1) 



where G is the Newtonian gravitational constant and c is the speed of light in vacuum. It acts on a 
test particle orbiting with a velocity v with the non-central acceleration [ |14| ] 

2 



ILT = VXfi g (1.2) 

c 



which induces secular precessions of the longitude of the ascending node Q. 



and the argument of pericentre co 
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'According to an interesting historical analysis recently performed by Pfister it would be more correct to speak 
about an Einstein-Thirring-Lense effect. 

2 It is another precessional effect predicted by general relativity ^ ^] which is not related to the mass-energy currents 
of the source. 
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Figure 1: Keplerian orbit. The longitude of the ascending node £2 is counted from a reference X direction 
in the equator of the central body, assumed as reference plane {X, Y}, to the line of the nodes which is the 
intersection of the orbital plane with the equatorial plane of the central body. It has mass M and proper 
angular momentum S. The argument of pericentre ft) is an angle in the orbital plane counted from the line 
of the nodes to the location of the pericentre, here marked with IT The time-dependent position of the 
moving test particle of mass m is given by the true anomaly /, counted anticlockwise from the pericentre's 
position. The inclination between the orbital and the equatorial planes is i. Thus, £2,0),/ can be viewed as the 
three (constant) Euler angles fixing the configuration of a rigid body, i.e. the orbit which in the unperturbed 
Keplerian case does change neither its shape nor its size, in the inertial {X, Y,Z} space. Courtesy by H.I.M. 
Lichtenegger, IWF, Graz. 




of the orbit of a test particle. In eq. ( |1.3| ) and eq. ( |1 .4| ) a and e are the semimajor axis and the 
eccentricity, respectively, of the test particle's orbit and i is its inclination to the central body's 
equator. The semimajor axis a fixes the size of the ellipse, while its shape is determined by the 
eccentricity < e < 1 ; an orbit with e = is a circle. The angles Q. and ft) establish the orientation 
of the orbit in the inertial space and in the orbital plane, respectively. O, ft) and i can be viewed as 
the three Euler angles which determine the orientation of a rigid body with respect to an inertial 
frame. In Figure [l] we illustrate the geometry of a Keplerian orbit [15]. 



In this Lecture we will critically discuss the following topics 

• Section ||. The realistic evaluation of the total accuracy in the test performed in recent years 
with the existing Earth 's artificial satellites LAGEOS and LAGEOS II [77[ [7<|/. 

LAGEOS was put into orbit in 1976, followed by its twin LAGEOS II in 1992; they are 
passive, spherical spacecraft entirely covered by retroreflectors which allow for their accurate 
tracking through laser pulses sent from Earth-based ground stations according to the Satellite 



Laser Ranging (SLR) technique [19]. They orbit at altitudes of about 6000 km (alageos = 
12270 km, #lageos ii = 12163 km) in nearly circular paths (^lageos = 0.0045, ^lageos ii = 
0.014) inclined by 110 deg and 52.65 deg, respectively, to the Earth's equator. The Lense- 
Thirring effect for their nodes amounts to about 30 milliarcseconds per year (mas yr _1 ) which 
correspond to about 1.7 m yr _1 in the cross-track direction 3 at the LAGEOS altitudes. 



3 A perturbing acceleration like Alt is customarily projected onto the radial r, transverse T and cross-track v direc- 
tions of an orthogonal frame comoving with the satellite [|l4t; it turns out that the Lense-Thirring node precession affects 
the cross-track component of the orbit according to Avlt ~ asiniAfJcr (eq. (A65), p. 6233 in Ref. [po|). 
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The idea of measuring the Lense-Thirring node rate with the just launched LAGEOS satel- 
lite, along with the other SLR targets orbiting at that time, was put forth by Cugusi and 
Proverbio ^2_ . Tests have started to be effectively performed later by using the LA- 



GEOS and LAGEOS II satellites [J23|] , according to a strategy by Ciufolini [g4J involving the 
use of a suitable linear combination of the nodes Q. of both satellites and the perigee (0 of 
LAGEOS II. This was done to reduce the impact of the most relevant source of systematic 
bias, i.e. the mismodelling in the even (£ = 2,4, 6 . . .) zonal (m = 0) harmonic coefficients 
Ji of the multipolar expansion of the Newtonian part of the terrestrial gravitational poten- 
tial due to the diurnal rotation (they induce secular precessions on the node and perigee of 
a terrestrial satellite much larger than the gravitomagnetic ones. The Ji coefficients cannot 
be theoretically computed but must be estimated by suitably processing long data sets from 
the dedicated satellites like CHAMP and GRACE; see Section |||): the three-elements com- 
bination used allowed for removing the uncertainties in J2 and J4. In Ref. [25] a « 20% test 
was reported by using the 4 EGM96 Earth gravity model [^] ; subsequent detailed analyses 
showed that such an evaluation of the total error budget was overly optimistic in view of the 
likely unreliable computation of the total bias due to the even zonals [|^, 28, 29]. An anal- 
ogous, huge underestimation turned out to hold also for the effect of the non-gravitational 



perturbations [ pQ| ] like the direct solar radiation pressure, the Earth's albedo, various subtle 
thermal effects depending on the physical properties of the satellites' surfaces and their ro- 



tational state [31, B2L 33, 34, 35 



,28], which the perigees of LAGEOS-like satellites 



are particularly sensitive to. As a result, the realistic total error budget in the test reported in 



[ 25 ] might be as large as 60 — 90% or (by considering EGM96 only) even more. 

The observable used 5 in Ref. pi] with the GRACE-only EIGEN-GRACE02S model [00] 



and in Ref. [18] with other more recent Earth gravity models was the following linear combi- 
nation 6 of the nodes of LAGEOS and LAGEOS II, explicitly computed in Ref. [ 42] following 



the approach put forth in Ref. [ J24Q 



f = a 



LAGEOS 



+ c l Cl LAGE0Sl \ 



where 



^LAGEOS 
" ^LAGEOS II 



COS ILAGEOS ( 1 - ^LAGEOS II 



COS/lagEOSII V i-^LAGEOS 



( ^LAGEOS II 
V ^LAGEOS 



7/2 



(1.5) 



(1.6) 



The coefficients Clj of the aliasing classical node precessions [ ]43[ ] £2 c iass = Y,i& £Ji induced 
by the even zonals have been analytically worked out up to I = 20 in, e.g., Ref. [p7[]; they 
yield c\ = 0.544. The Lense-Thirring signature of eq. (1.5) amounts to 47.8 mas yr _1 . 



The combination eq. (1.5) allows, by construction, to remove the aliasing effects due to 



4 Contrary to the subsequent models based on the dedicated satellites CHAMP (http://www- 
app2.gfz-potsdam.de/pbl/op/champ/index_CHAMP.html) and GRACE (http://www-app2.gfz- 

potsdam.de/pbl/op/grace/index_GRACE.html), EGM96 relies upon multidecadal tracking of SLR data of a con- 
stellation of geodetic satellites including LAGEOS and LAGEOS II as well; thus the possibility of a sort of a — priori 
'imprinting' of the Lense-Thirring effect itself, not solved-for in EGM96, cannot be neglected. 

5 The first test was, in fact, performed with the CHAMP-only based model EIGEN-2 by Lucchesi in Ref. 

6 See also Refs. @ || |§]. 
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the static and time-varying parts of the first even zonal J2. The nominal (i.e. computed 
with the estimated values of Je, £ = 4,6...) bias due to the remaining higher degree even 
zonals would amount to about 10 5 mas yr _1 ; the need of a careful and reliable modeling of 
such an important source of systematic bias is, thus, quite apparent. Conversely, the nodes 
of the LAGEOS-type spacecraft are directly affected by the non-gravitational accelerations 
1% level of the Lense-Thirring effect [|33|, 34, 35, 36, 37]. For a comprehensive, 



at a 



up-to-date overview of the numerous and subtle issues concerning the measurement of the 



Lense-Thirring effect see, e.g., Ref. [44]. 

Section [§. Another approach which could be followed in extracting the Lense-Thirring effect 
from the data of the LAGEOS-type satellites. 

Section |j. The possibility that the LARES mission, recently approved by the Italian Space 
Agency (ASI), will be able to measure the Lense-Thirring node precession with an accuracy 
of the order of\%. 



In Refs. [|45j, |46J] it was proposed to measure the Lense-Thirring precession of the nodes D. 
of a pair of counter-orbiting spacecraft to be launched in terrestrial polar orbits and endowed 
with drag-free apparatus. A somewhat equivalent, cheaper version of such an idea was put 



forth in 1986 by Ciufolini in Ref. [ |47[ ] who proposed to launch a passive, geodetic satellite 
in an orbit identical to that of LAGEOS apart from the orbital planes which should have 
been displaced by 180 deg apart. The measurable quantity was, in the case of the proposal of 
Ref. p7|], the sum of the nodes of LAGEOS and of the new spacecraft, later named LAGEOS 
III, LARES, WEBER-SAT, in order to cancel to a high level of accuracy the corrupting 
effect of the multipoles of the Newtonian part of the terrestrial gravitational potential which 
represent the major source of systematic error (see Section |2|). Although extensively studied 

it 



by various groups [ |48| , |49j ], such an idea was not implemented for many years. In Ref. 
was proposed to include also the data from LAGEOS II by using a different observable. Such 
an approach was proven in Ref. to be potentially useful in making the constraints on the 
orbital configuration of the new SLR satellite less stringent than it was originally required 
in view of the recent improvements in our knowledge of the classical part of the terrestrial 
gravitational potential due to the dedicated CHAMP and, especially, GRACE missions. 

Since reaching high altitudes and minimizing the unavoidable orbital injection errors is ex- 
pensive, it was explored the possibility of discarding LAGEOS and LAGEOS II using a 



low-altitude, nearly polar orbit for LARES Q52| , p3|], but in Refs. [|5_4J, |55|] it was proven that 



such alternative approaches are not feasible. It was also suggested that LARES would be 



able to probe alternative theories of gravity [56], but also in this case it turned out to be 



impossible []57|, |58|]. 

The stalemate came to an end when about one year ago ASI made the following official 
announcement 7 : "On February 8, the ASI board approved funding for the LARES mis- 
sion, that will be launched with VEGAs maiden flight before the end of 2008. LARES is 
a passive satellite with laser mirrors, and will be used to measure the Lense-Thirring ef- 
fect." The italian version of the announcement yields some more information specifying 



7 It was posted on the currently expired webpage http://www.asi.it/SiteEN/MotorSearchFullText.aspx?keyw=LARES. 
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that LARES, designed in collaboration with National Institute of Nuclear Physics (INFN), 
is currently under construction by Carlo Gavazzi Space SpA; its Principal Investigator (PI) 
is I. Ciufolini and its scientific goal is to measure at a 1% level the Lense-Thirring effect 
in the gravitational field of the Earth. Concerning the orbital configuration of LARES, in 
one of the latest communication to INFN, Rome, 30 January 2008, Ciufolini [ [59^ ] writes 
that LARES will be launched with a semimajor axis of approximately 7600 km and an 
inclination between 60 and 80 deg. More precise information can be retrieved in Sec- 
tion 5.1, pag 9 of the document Educational Payload on the Vega Maiden Flight Call For 
CubeSat Proposals, European Space Agency, Issue 111 February 2008, downloadable at 
http://esamultimedia.esa.int/docs/LEX-EC/CubeSat%20CFP%20issue%201.pdf. It is writ- 
ten there that LARES will be launched into a circular orbit with altitude h = 1200 km, 
corresponding to a semimajor axis czlares = 7578 km, and inclination i = 71 deg to the 
Earth's equator. Latest information 8 point towards a launch at the end of 2009 with a 
VEGA rocket in a circular orbit inclined by 71 deg to the Earth's equator at an altitude 
of 9 1450 km corresponding to a semimajor axis of gilr = 7828 km. More or less the same 
has been reported by Ciufolini to the INFN in Villa Mondragone, 3 October 2008 [|50|], 
and to the delegates of the 16th International Workshop on Laser Ranging, 13-17 October 
2008 [^l|]. However, according to more recent information from the official ASI website 
(http://www.asi.it/en/activity/cosmology/lares), it seems that there is still a broad uncertainty 
in the final orbital configuration of LARES which should, now, correspond to a as 7778 km 
and 60 deg < i < 86 deg. 

Section || The possibility of measuring the Lense-Thirring effect of the Sun with the inner 
planets of the Solar System. 

Recent determinations of the Sun's proper angular momentum 

S = (19O.O±1.5)xlO 39 kgm 2 s- 1 (1.7) 

from helioseismology [p2,|63], accurate to 0.8%, yields a value about one order of magnitude 
smaller than that obtained by assuming an homogeneous and uniformly rotating Sun, as done 
in the pionieristic work by De Sitter [64], and also in Refs. [|l4], ^1 who concluded that, at 



their time, it was not possible to measure the solar Lense-Thirring effect. Now, instead, 
the situation seems much more promising, in spite of the reduced magnitude of the solar 
gravitomagnetic field with respect to the earlier predictions. As in the Schwarzschild metric 
the characteristic length to be compared with the accuracy in determining particles 'orbits is 
the Schwarzschild radius R g = 2GM/c = 3 km for the Sun, the typical solar gravitomagnetic 
length is 

l ? = T7^= 3l9m - ( Lg ) 

g MqC 



8 See on the WEB http://www.esa.int/esapub/bulletin/bulletinl35/bull35f_bianchi.pdf. 

9 I thank Dr. D. Barbagallo (ESRIN) for having kindly provided me with the latest details of the orbital configuration 



of LARES. 
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Table 1: First line: uncertainties, in m, in the average heliocentric distances of the inner planets obtained by 
propagating the formal errors in a and e according to Table 3 of Ref. [|65p; the EPM2006 ephemerides were 
used by E.V. Pitjeva [|65||. Second line: maximum differences, in m, between the EPM2006 and the DE414 
[ |56| l ephemerides for thr inner planets in the time interval 1 960-2020 according to Table 5 of Ref. [ |&7| ] . They 
have to be compared to the characteristic gravitomagnetic length of the Sun Z® = 319 m. 



Type of orbit uncertainty Mercury Venus Earth Mars 

8 (r) (EPM2006) 38 3 1 2 

EPM2006-DE414 256 131 17.2 78.7 



( r )=a[l + e -\, (1.9) 



The present-day accuracy in knowing, e.g., the inner planet's mean radius 

2 

is shown in Table ffl. Such values have been obtained by linearly propagating the formal, sta- 
tistical errors in a and e according to Table 3 of Ref. [ |6"5j ] ; it is clear that, even by re-scaling 
them by a a factor of, say, 2 — 5, the gravitomagnetic effects due to the Sun's rotation fall 
in the measurability domain. Another possible way to evaluate the present-day uncertainty 
in the planetary orbital motions consists of looking at different ephemerides of comparable 
accuracy. In Table |l] we do that for the EPM2006 [|5]| and the DE414 ||6|] ephemerides; 



although, larger than 8 (r), the maximum differences between such ephemerides are smaller 
than the solar gravitomagnetic length Z®. Since much ranging data to Venus Express space- 
craft 10 will be collected and processed, and in view of the ongoing Messenger [ ]68| ] and the 
future Bepi-Colombo [ ]69| , |68| ] missions to Mercury and of the developments in the Planetary 
Laser Ranging (PLR) technique J7C|], it seems that the planetary orbit accuracy will be further 
improved. It is remarkable to note that the currently available estimate of S Q from helioseis- 
mology is accurate enough to allow, in principle, a genuine Lense-Thirring test. Moreover, it 
was determined in a relativity-free fashion from astrophysical techniques which do not rely 
on the dynamics of planets in the gravitational field of the Sun. Thus, there is no any a priori 
"memory" effect of general relativity itself in the adopted value of S Q . 

• Section ^. The test with the Mars Global Surveyor probe in the gravitational field of Mars. 
Since the angular momentum of Mars can be evaluated as 

S M = (1.92±0.01) x 10 32 kgm 2 s" 1 (1.10) 

from the latest spacecraft-based determinations of the areophysical parameters quoted in 
Ref. 0], it turns out 

|M = ^L = LOm , (U1) 
8 M M c 

Such a value has to be compared with the present-day accuracy in determining the orbit of a 
spacecraft like, e.g., Mars Global Surveyor (MGS); it can be evaluated it is about 0.15 m 



10 



See on the WEB http://www.esa.int/SPECIALS/Venus_Express/. 
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in the radial direction, not affected by the gravitomagnetic force itself. Thus, it makes sense 
to investigate the possibility of measuring the Lense-Thirring effect in the gravitational field 
of Mars as well. In fact, the time series of the Root Mean Square (RMS) orbit overlap differ- 
ences 0] of the out-of-plane part v of the orbit of the nearly polar MGS probe (a = 3796 
km, i = 92.86 deg, e = 0.0085) over a time span AP « 5 yr, processed without modelling 
gravitomagnetism at all, has been successfully interpreted in terms of the Lense-Thirring ef- 



fect by Iorio in Refs. [72, 73]. However, criticisms have been raised by Krogh in Ref. [74]; 



Iorio has replied in Ref. [173]. In Section || we will discuss them in details. 
2. The systematic error of gravitational origin in the LAGEOS-LAGEOS II test 



The realistic evaluation of the total error budget of the LAGEOS-LAGEOS II node test [16] 



raised a lively debate [76, 17, 51, 77, 78, [7^, 8C, 81], mainly focussed on the impact of the static 
and time-varying parts of the Newtonian component of the Earth's gravitational potential through 
the secular precessions induced on a satellite's node. 

In the real world the path of a probe is not only affected by the relativistic gravitomagentic field 
but also by a huge number of other competing classical orbital perturbations of gravitational and 
non-gravitational origin. The most insidious disturbances are those induced by the static part of the 
Newtonian component of the multipolar expansion in spherical harmonics 11 j£,£ = 2,4,6, ... of the 



gravitational potential of the central rotating mass [ ]43[ ] : they affect the node with effects having the 
same signature of the relativistic signal of interest, i.e. linear trends which are orders of magnitude 
larger and cannot be removed from the time series of data without affecting the Lense-Thirring 
pattern itself as well. The only thing that can be done is to model such a corrupting effect as most 
accurately as possible and assessing the impact of the residual mismodelling on the measurement of 
the frame-dragging effect. The secular precessions induced by the even zonals of the geopotential 
can be written as 

Qgecpot = £^, (21) 

where the coefficients £l.g,£ = 2,4,6, ... depend on the parameters of the Earth (GM and the equa- 
torial radius R) and on the semimajor axis a, the inclination i and the eccentricity e of the satellite. 
For example, for I = 2 we have 

3 / R \ 2 cos i 

n = ^GMja? is the Keplerian mean motion. They have been analytically computed up to i = 20 



in, e.g., [27]. Their mismodelling can be written as 



8& e °v ot < J^\Clj\SJ e , (2.3) 

e=2 



where SJt represents our uncertainty in the knowledge of the even zonals J% 



The relation among the even zonals J( and the normalized gravity coefficients Qrj which are customarily deter- 
mined in the Earth's gravity models, is Jg = —y*2E+ 1 Qq. 
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A common feature of all the competing evaluations so far published is that the systematic 
bias due to the static component of the geopotential was always calculated by using the released 
(more or less accurately calibrated) sigmas Oj ( of one Earth gravity model solution at a time for the 
uncertainties 8J^. Thus, it was said that the model X yields a x% error, the model Y yields a y% 
error, and so on. 

Since a unstable calibration of the formal, statistical uncertainties in the estimated zonals of 
the covariance matrix of a global solution is always a difficult task to be implemented in a reliable 
way, a much more realistic and conservative approach consists, instead, of taking the difference 12 



A7/ = |4(X)-J*(Y)|, 1 = 2,4,6, 



(2.4) 



of the estimated even zonals for different pairs of Earth gravity field solutions as representative of 
the real uncertainty 8 J? in the zonals [33]. In Tables-Table [10| we present our results for the most 
recent GRACE-based models released so far by different institutions and retrievable on the Internet 
at 13 http://icgem.gfz-potsdam.de/ICGEM/ICGEM.html. The models used are EIGEN-GRACE02S 
[@ from GFZ (Potsdam, Germany), GGM02S [Q and GGM03S [||] from CSR (Austin, Texas), 



ITG-Grace02s [|86J and ITG-Grace03s [g7j] from IGG (Bonn, Germany), JEM01-RL03B from JPL 



(NASA, USA) and AIUB-GRACE01S []88|] from AIUB (Switzerland). Note that this approach was 
explicitly followed also by Ciufolini in Ref. [||] with the JGM3 and GEMT-2 models. In Table 
^|-Table |lO| we quote both the sum Y,e=4fi °f tne absolute values of the individual mismodelled 
terms 

ft = \Cl L £ AGEOS + d^ AGEOS 11 1 AJg (2.5) 



(SAV), and the square root of the sum of their squares y X^Lt/f (RSS); in both cases we normal- 
ized them to the combined Lense-Thirring total precession of 47.8 mas yr _1 . 

The systematic bias evaluated with a more realistic approach is about 3 to 4 times larger than 
one can obtain by only using this or that particular model. The scatter is still quite large and far 
from the 5 — 10% claimed in Ref. [16]. In particular, it appears that J4, and to a lesser extent 



78, which are just the most relevant zonals for us because of their impact on the combination of 



eq. (|L5|), are the most uncertain ones, with discrepancies AJg between different models, in general, 
larger than the sum of their sigmas Oj n calibrated or not. 

Another way to evaluate the uncertainty in the LAGEOS-LAGEOS II node test may consist of 
computing the nominal values of the total combined precessions for different models and compar- 
ing them, i.e. by taking 



^LAGEOS + Cl ^LAGEOS II) ^ (x) _ j^y)] 



/=4 



(2.6) 



The results are shown in Table 1 1 



A different approach that could be followed to take into account the scatter among the various 
solutions consists in computing mean and standard deviation of the entire set of values of the 



12 See Fig.5 of [|^] for a comparison of the estimated C40 in different models. 

13 I thank M Watkins (JPL) for having provided me with the even zonals and their sigmas of the JEM01-RL03B 
model. 
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Table 



LAGEOS 



Impact of the 

■ Cl Qh mmsn \AJ e , £ = 



mismodelling in the even zonal 
4,..., 20, in mas yr -1 . Recall that J( 



harmonics on fy = 
= -VW+T Cro; for the 
between the model 



"WO 



uncertainty in the even zonals we have taken here the difference AQo 

X=GGM02S [0 and the model Y=ITG-Grace02s[|6). GGM02S is based on 363 days of GRACE-only 
data (GPS and intersatellite tracking, neither constraints nor regularization applied) spread between April 4, 
2002 and Dec 31, 2003. The a are formal for both models. AQo are always larger than the linearly added 
sigmas, apart from £ — 12 and £ = 18. Values of smaller than 0.1 mas yr~' have not been quoted. The 
Lense-Thirring precession of the combination of eq. (L5) amounts to 47.8 mas yr . The percent bias 8ji 
have been computed by normalizing the linear sum of fg, £ = 4, ... ,20 (S AV) and the square root of the sum 
of ff,l — 4, ... ,20 to the Lense-Thirring combined precessions. 



1 


AQo (GGM02S-ITG-Grace02s) 


Cx + CTy 


ft (mas yr ! ) 


4 


1.9 x lO" 11 


8.7 x 10~ 12 


7.2 


6 


2.1 x lO" 11 


4.6 x 10~ 12 


4.6 


8 


5.7 x 10~ 12 


2.8 x 10~ 12 


0.2 


10 


4.5 x 10~ 12 


2.0 x 10~ 12 




12 


1.5 x 10~ 12 


1.8 x 10~ 12 




14 


6.6 x 10~ 12 


1.6 x 10~ 12 




16 


2.9 x 10~ 12 


1.6 x 10~ 12 




18 


1.4 x 10~ 12 


1.6 x 10~ 12 




20 


2.0 x 10~ 12 


1.6 x 10~ 12 






8 =25% (SAV) 


8pL = 18% (RSS) 





Table 3: Bias due to the mismodelling in the even zonals of the models X=ITG-Grace03s [p7|], based on 
GRACE-only accumulated normal equations from data out of September 2002-April 2007 (neither apriori 
information nor regularization used), and Y=GGM02S fl84}]. The a for both models are formal. AQo are 
always larger than the linearly added sigmas, apart from £=12 and £ = 18. 



1 


AQo (ITG-Grace03s-GGM02S) 


Cx + cty 


ft (masyr ! ) 


4 


2.58 x lO" 11 


8.6 x 10~ 12 


9.6 


6 


1.39 x lO" 11 


4.7 x 10" 12 


3.1 


8 


5.6 x 10~ 12 


2.9 x 10~ 12 


0.2 


10 


1.03 x 10- 11 


2x 10- 12 




12 


7 x 10~ 13 


1.8 x 10" 12 




14 


7.3 x 10~ 12 


1.6 x 10~ 12 




16 


2.6 x 10~ 12 


1.6 x 10- 12 




18 


8 x 10~ 13 


1.6 x 10~ 12 




20 


2.4 x 10~ 12 


1.6 x 10" 12 






b\i = 27% (SAV) 


Sn =21% (RSS) 
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Table 4: Bias due to the mismodelling in the even zonals of the models X = GGM02S @ and Y = GGM03S 
[ p5| | retrieved from data spanning January 2003 to December 2006. The <T for GGM03S are calibrated. AQo 
are larger than the linearly added sigmas for £ = 4,6. (The other zonals are of no concern) 







Ox + Oy 


ft ( mas y r ! ) 


4 


1.87 x 10~ n 


1.25 x lO" 11 


6.9 


6 


1.96 x 10~ n 


6.7 x 10~ 12 


4.2 


8 


3.8 x 10~ 12 


4.3 x 10~ 12 


0.1 


10 


8.9 x 10~ 12 


2.8 x 10~ 12 


0.1 


12 


6x 10~ 13 


2.4 x 10~ 12 




14 


6.6 x 10~ 12 


2.1 x 10~ 12 




16 


2.1 x 10~ 12 


2.0 x 10~ 12 




18 


1.8 x 10~ 12 


2.0 x 10~ 12 




20 


2.2 x 10~ 12 


1.9 x 10~ 12 






8 pL = 24% (SAV) 


8}X = 17% (RSS) 





Table 5: Bias due to the mismodelling in the even zonals of the models X = EIGEN-GRACE02S and 
Y = GGM03S p5|]. The a for both models are calibrated. AQo are always larger than the linearly added 
sigmas apart from £ = 14, 18. 



i 


AQo (EIGEN-GRACE02S -GGM03S ) 


Ox + Oy 


ft (masyr *) 


4 


2.00 x 10- 11 


8.1 x 10~ 12 


7.4 


6 


2.92 x lO" 11 


4.3 x 10~ 12 


6.3 


8 


1.05 x lO" 11 


3.0 x 10~ 12 


0.4 


10 


7.8 x 10~ 12 


2.9 x 10~ 12 


0.1 


12 


3.9 x 10~ 12 


1.8 x 10~ 12 




14 


5 x 10~ 13 


1.7 x 10~ 12 




16 


1.7 x 10~ 12 


1.4 x 10~ 12 




18 


2x 10~ 13 


1.4 x 10~ 12 




20 


2.5 x 10~ 12 


1.4 x 10~ 12 






8\i = 30% (SAV) 


8pL = 20% (RSS) 





even zonals for the models considered so far, degree by degree, and taking the standard deviations 
as representative of the uncertainties 8Je,£ = 4,6,8,.... It yields 8jJ. = 15%, in agreement with 



Ref. [18] 



It must be recalled that also the further bias due to the cross-coupling between Ji and the orbit 



inclination, evaluated to be about 9% in Ref. [79], must be added 
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Table 6: Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B, based on 49 
months of GRACE-only data, and Y = GGM03S The a for GGM03S are calibrated. AC f0 are always 
larger than the linearly added sigmas apart from I = 16. 



t 


AQo ( JEMO 1 -RL03B -GGM03 S) 


Ox + CTy 


ft (mas yr 1 ) 


4 


i.y / x iu 


A O *y 1 f\— 12 

4. J X 1U 


f.J 


6 


2.7 x 10~ 12 


2.3 x 10~ 12 


0.6 


8 


1.7 x 10~ 12 


1.6 x 10~ 12 




10 


2.3 x 10~ 12 


8 x 10~ 13 




12 


7 x 10~ 13 


7 x 10~ 13 




14 


1.0 x 10~ 12 


6x 10~ 13 




16 


2x 10~ 13 


5 x 10~ 13 




18 


7 x 10~ 13 


5 x 10~ 13 




20 


5 x 10~ 13 


4x 10~ 13 






8/U = 17% (SAV) 


8n = 15% (RSS) 





Table 7: Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B and Y = ITG- 
Grace03s [p7|]. The a for ITG-Grace03s are formal. AQo are always larger than the linearly added sigmas. 



1 


AQo (JEM01-RL03B-ITG-Grace03s) 


Cx + Cy 


ft (mas yr ! ) 


4 


2.68 x 10~ n 


4x 10~ 13 


9.9 


6 


3.0 x 10~ 12 


2x 10~ 13 


0.6 


8 


3.4 x 10~ 12 


1 x 10~ 13 


0.1 


10 


3.6 x 10~ 12 


1 x 10~ 13 




12 


6x 10~ 13 


9 x 10~ 14 




14 


1.7 x 10~ 12 


9x 10~ 14 




16 


4x 10~ 13 


8 x 10~ 14 




18 


4x 10~ 13 


8 x 10~ 14 




20 


7 x 10~ 13 


8 x 10~ 14 






5jU = 22% (SAV) 


<5jU = 10% (RSS) 





3. Possible alternative approaches to extract the Lense-Thirring signature from the 
data 



The technique adopted so far in Ref. [16] and Ref. [hq] to extract the gravitomagentic signal 



from the LAGEOS and LAGEOS II data is described in detail in, e.g., Refs. [89, 



In both the approaches the action of the gravitomagnetic field is not included in the dynamical 
force models used to fit the satellites' data. In the data reduction process no dedicated grav- 
itomagnetic parameter is estimated, contrary to, e.g., station coordinates, state vector, satel- 
lites' drag and radiation coefficients Cp and Cr, respectively, etc.; the Lense-Thirring pattern 
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Table 8: Aliasing effect of the mismodelling in the even zonal harmonics estimated in the X=ITG-Grace03s 
10 and the Y=EIGEN-GRACE02S |@ models. The covariance matrix a for ITG-Grace03s are formal, 
while the ones of EIGEN-GRACE02S are calibrated. AQo are larger than the linearly added sigmas for 
£ = 4,.. .,20, apart from 1= 18. 



t 


AQo (ITG-Grace03s-EIGEN-GRACE02S) 


0x + <Jy 


ft (mas yr l ) 


4 


2.72 x lO" 11 


3.9 x 10~ 12 


10.1 


6 


2.35 x 10~ n 


2.0x 10~ 12 


5.1 


8 


1.23 x 10- 11 


1.5 x 10~ 12 


0.4 


10 


9.2 x 10~ 12 


2.1 x 10~ 12 


0.1 


12 


4.1 x 10~ 12 


1.2 x 10~ 12 




14 


5.8 x 10~ 12 


1.2x 10~ 12 




16 


3.4 x 10~ 12 


9x 10~ 13 




18 


5 x 10~ 13 


1.0 x 10~ 12 




20 


1.8 x 10~ 12 


1.1 x 10~ 12 






8 n = 37% (SAV) 


SjU = 24% (RSS) 





Table 9: Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B, based on 49 
months of GRACE-only data, and Y = AIUB-GRACE01S [[|§. The latter one was obtained from GPS 
satellite-to-satellite tracking data and K-band range-rate data out of the period January 2003 to December 
2003 using the Celestial Mechanics Approach. No accelerometer data, no de-aliasing products, and no 
regularisation was applied. The a for AIUB-GRACE01S are formal. AQo are always larger than the linearly 
added sigmas. 



£ 


AQo ( JEM0 1 -RL03B - AIUB -GR ACE0 1 S ) 


Gx + CTy 


/Hmasyr l ) 


4 


2.95 x 10~ n 


2.1 x 10~ 12 


11 


6 


3.5 x 10~ 12 


1.3 x 10~ 12 


0.8 


8 


2.14 x 10- 11 


5 x 10~ 13 


0.7 


10 


4.8 x 10~ 12 


5 x 10~ 13 




12 


4.2 x 10~ 12 


5 x 10~ 13 




14 


3.6 x 10~ 12 


5 x 10~ 13 




16 


8 x 10~ 13 


5 x 10~ 13 




18 


7 x 10~ 13 


5 x 10~ 13 




20 


1.0 x 10~ 12 


5 x 10~ 13 





8pL = 26% (SAV) SjU = 23% (RSS) 



is retrieved with a sort of post-post-fit analysis in which the time series of the computed 
"residuals" of the nodes with the difference between the orbital elements of consecutive arcs, 



The expression "residuals of the nodes" is used, strictly speaking, in an improper sense because the Keplerian 
orbital elements are not directly measured. 
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Table 10: Bias due to the mismodelling in the even zonals of the models X = EIGEN-GRACE02S |@ and 
Y = AIUB-GRACE01S The a for AIUB-GRACE01S are formal, while those of EIGEN-GRACE02S 
are calibrated. AQy are larger than the linearly added sigmas for £ = 4,6,8, 16. 



1 


AQ (EIGEN-GRACE02S — AIUB -GR ACEO 1 S ) 


ox + oy 


fp (mas yr l ) 


4 


2.98 x 10~ n 


6.0 x 10 12 


11.1 


6 


2.29 x lO" 11 


3.3 x 10- 12 


5.0 


8 


1.26 x 10- 11 


1.9 x 10~ 12 


0.4 


10 


6x 10~ 13 


2.5 x 10~ 12 




12 


5 x 10~ 13 


1.6 x 10~ 12 




14 


5 x 10~ 13 


1.6 x 10~ 12 




16 


2.9 x 10~ 12 


1.4 x 10~ 12 




18 


6x 10~ 13 


1.4 x 10~ 12 




20 


2x 10~ 13 


1.5 x 10- 12 






Sfi = 34% (SAV) 


SjU = 25% (RSS) 





combined with eq. (1.5), is fitted with a straight line 



15 



In regard to possible other approaches which could be followed, it would be useful to, e.g., es- 
timate (in the least square sense), among other solve-for parameters, purely phenomenologi- 
cal corrections A& to the LAGEOS/LAGEOS II node precessions as well, without modelling 
the Lense-Thirring itself, so that it will be, in principle, contained in Af2, and combine them 



according to eq. (|L5|). Something similar has been done-although for different scopes-for the 
perihelia of the inner planets of the Solar System and the periastron of the pulsars [pl|]. 
To be more definite, various solutions with a complete suite of dynamical models, apart from 
the gravitomagnetic force itself, should be produced in which one inserts a further solve-for 
parameter, i.e. a correction A& to the standard Newtonian modelled precessions. One could 
see how the outcome varies by changing the data sets and/or the parameters to be solved for. 
Maybe it could be done for each arc, so to have a collection of such node extra-rates. Such a 
strategy would be much more model-independent. 

• Another previously suggested [^] way to tackle the problem consists of looking at a Lense- 
Thirring-dedicated parameter to be estimated along with all the zonals in a new global solu- 
tion for the gravity field incorporating the gravitomagnetic component as well; instead, in all 
the so far produced global gravity solutions no relativistic parameter(s) have been included 
in the set of the estimated ones. 

A first, tentative step towards the implementation of the strategy of the first point mentioned 
above with the LAGEOS satellites in terms of the PPN parameter y has been recently taken by 
Combrinck in Ref. I^l. 



15 Plus some harmonic signals as well to remove the time-dependent tidal perturbations. 
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Table 11: Systematic uncertainty 8ji in the LAGEOS-LAGEOS II test evaluated by taking the absolute 
value of the difference between the nominal values of the total combined node precessions due to the even 
zonals for different models X and Y, i.e. |ns eo P ot (X) - n§ eo P ot (Y) I . 



Models compared 8 /I 

AIUB -GR ACEO IS— JEMO 1 -RL03B 20% 

AIUB -GR ACEO IS- GGM02S 27% 

AIUB -GR ACEO IS- GGM03S 3% 

AIUB-GRACE01S-ITG-Grace02 2% 

AIUB-GRACE01S-ITG-Grace03 0.1% 
AIUB -GR ACEO 1 S — EIGEN-GRACE02S 33% 

JEM01 -RL03B — GGM02S 7% 

JEM01-RL03B— GGM03S 17% 

JEM01-RL03B-ITG-Grace02 18% 

JEM01-RL03B-ITG-Grace03s 20% 

JEM01 -RL03B — EIGEN-GRACE02S 13% 

GGM02S — GGM03S 24% 

GGM02S-ITG-Grace02 25% 

GGM02S-ITG-Grace03s 27% 

GGM02S — EIGEN-GRACE02S 6% 

GGM03S-ITG-Grace02 1% 

GGM03S-ITG-Grace03s 3% 

GGM03S — EIGEN-GRACE02S 30% 

ITG-Grace02-ITG-Grace03s 2% 

ITG-Grace02 — EIGEN-GRACE02S 3 1 % 

ITG-Grace03s— EIGEN-GRACE02S 33% 



4. On the LARES mission 

The combination that should be used for measuring the Lense-Thirring effect with LAGEOS, 
LAGEOS II and LARES is [0 

^LAGEOS + ^LAGEOS II + ^LARES _ (4 , } 



The coefficients k[ and &2 entering eq. (4.1) are defined as 



ALARES ALAGEOS _ ALAGEOS ALARES 

/, 11 4 11 2 4 A 'ICQA 

— fiLAGEOS II ALARES _qLARES£jLAGEOS II — U.-3JOU, 



ALAGEOS ALAGEOS II ALAGEOS II ALAGEOS 

?,„ "2 ".4 i ±2 ".4 n C\1K 1 

K 2 — aLaGeOs UaLArEs aLareS aLaGLos H — U.U/31. 



(4.2) 



The combination of eq. (4.1) cancels out, by construction, the impact of the first two even zonals; 
we have used cjlr = 7828 km, /lr = 71.5 deg. The total Lense-Thirring effect, according to eq. 
( p~T| ) and eq. (4.2), amounts to 50.8 mas yr _1 . 
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4.1 A conservative evaluation of the impact of the geopotential on the LARES mission 

The systematic error due to the uncancelled even zonals J6,J&,... can be conservatively evalu- 
ated as 

SH < £ | ^ AGEOS + ^ AGEOS 11 + k 2 n L £ ARES | 8J ( (4.3) 

1=6 

Of crucial importance is how to assess 8J{. By proceeding as in Section |2| and by using the 
same models up to degree £ = 60 because of the lower altitude of LARES with respect to LAGEOS 



and LAGEOS II which brings into play more even zonals, we have the results presented in Table [12. 
They have been obtained with the standard and widely used Kaula approach [ |43| ] in the following 
way. We, first, calibrated our numerical calculation with the analytical ones performed with the 
explicit expressions for Cl g worked out up to £ = 20 in Ref. [|27|]; then, after having obtained 
identical results, we confidently extended our numerical calculation to higher degrees by means of 
two different softwares 



Table 12: Systematic percent uncertainty 8ji in the combined Lense-Thirring effect with LAGEOS, LA- 
GEOS II and LARES according to eq. (4.3) and 5J( = AJ( up to degree £ — 60 for the global Earth's 
gravity solutions considered here; the approach by Kaula j43| ] has been followed. For LARES we adopted 
«lr = 7828 km, /lr = 71.5 deg, <?lr = 0.0. 



Models compared (8J e = AJ e ) Sfi (SAV) 8pL (RSS) 



AIUB -GRACE0 1 S - JEM0 1 -RL03B 


23% 


16% 


AIUB -GRACE0 1 S -GGM02S 


16% 


8% 


AIUB -GRACE0 1 S -GGM03 S 


22% 


13% 


AIUB -GRACE0 1 S - ITG-Grace02 


24% 


15% 


AIUB -GRACE0 1 S - LTG-Grace03 


22% 


14% 


AIUB -GRACE01 S — EIGEN-GRACE02S 


14% 


7% 


JEM0 1 -RL03B - GGM02S 


14% 


9% 


JEM0 1 -RL03B -GGM03S 


5% 


3% 


JEM0 1 -RL03B - ITG-Grace02 


4% 


2% 


JEM01-RL03B-ITG-Grace03s 


5% 


2% 


JEM0 1 -RL03B — EIGEN-GRACE02S 


26% 


15% 


GGM02S— GGM03S 


13% 


7% 


GGM02S - ITG-Grace02 


16% 


8% 


GGM02S-ITG-Grace03s 


14% 


7% 


GGM02S -EIGEN-GRACE02S 


14% 


7% 


GGM03S-ITG-Grace02 


3% 


2% 


GGM03S-ITG-Grace03s 


2% 


0.5% 


GGM03S -EIGEN-GRACE02S 


24% 


13% 


ITG-Grace02-ITG-Grace03s 


3% 


2% 


ITG-Grace02— EIGEN-GRACE02S 


25% 


14% 


ITG-Grace03s — EIGEN-GRACE02S 


24% 


13% 
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It must be stressed that they may be still optimistic: indeed, computations for I > 60 become 
unreliable because of numerical instability of the results. 



In Table |3j we repeat the calculation by using for dJe the covariance matrix sigmas Oj e ; also 

60. 



in this case we use the approach by Kaula [ ]43[ ] up to degree 



Table 13: Systematic percent uncertainty 8/1 in the combined Lense-Thirring effect with LAGEOS, LA- 
GEOS II and LARES according to eq. (4.3) and 8Jf = 0>, up to degree £ = 60 for the global Earth's 
gravity solutions considered here; the approach by Kaula J43| ] has been followed. For LARES we adopted 
a LR = 7828 km, / LR = 71.5 deg, e LR = 0.0. 



Model (8Ji = a t ) <SjU (SAV) 8pL (RSS) 



AIUB-GRACE01S (formal) 


11% 


9% 


JEM01-RL03B 


1% 


0.9% 


GGM03S (calibrated) 


5% 


4% 


GGM02S (formal) 


20% 


15% 


ITG-Grace03s (formal) 


0.3% 


0.2% 


ITG-Grace02s (formal) 


0.4% 


0.2% 


EIGEN-GRACE02S (calibrated) 


21% 


17% 



If, instead, one assumes 8J( = se, £ = 2,4,6, ... i.e., the standard deviations of the sets of all 
the best estimates of J( for the models considered here the systematic bias, up to £ = 60, amounts 
to 12% (SAV) and 6% (RSS). Again, also this result may turn out to be optimistic for the same 
reasons as before; the analysis by Iorio of Ref. [ p5j ] seems to confirm such concerns. 

It must be pointed out that the evaluations presented here rely upon calculations of the coeffi- 



cients Q.£ performed with the well known standard approach by Kaula [43]; it would be important 



to try to follow also different computational strategies in order to test them. 

4.2 The impact of some non-gravitational perturbations 

It is worthwhile noting that also the impact of the subtle non-gravitational perturbations will 
be different with respect to the original proposal because LARES will fly in a lower orbit and 
its thermal behavior will probably be different with respect to LAGEOS and LAGEOS II. The 
reduction of the impact of the thermal accelerations, like the Yarkovsky-Schach effects, should 



have been reached with two concentric spheres. However, as explained in Ref. [|96|], this solution 
will increase the floating potential of LARES because of the much higher electrical resistivity and, 
thus, the perturbative effects produced by the charged particle drag. Moreover, the atmospheric 
drag will increase also because of the lower orbit of the satellite, both in its neutral and charged 
components. Indeed, although it does not affect directly the node O, it induces a secular decrease 



of the inclination i of a LAGEOS -like satellite [30] which translates into a further bias for the node 
itself according to 

„ ■ 3 / R\ 2 sin/ 7? „ 

M tag = ^(-J j^Si, (4.4) 
in which Si accounts not only for the measurement errors in the inclination, but also for any un- 



modelled/mismodelled dynamical effect on it. According to Ref. [97], the secular decrease for 
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LARES would amount to 

dt/ 



=y -0.6 mas yr 1 (4.5) 

LR 



yielding a systematic uncertainty in the Lense-Thirring signal of eq. (1.1) of about 3 — 9% yr 



An analogous indirect node effect via the inclination could be induced by the thermal Yarkovski- 



Rubincam force as well [97]. Also the Earth's albedo, with its anisotropic components, may have 
a non-negligible effect. 

Let us point out the following issue as well. At present, it is not yet clear how the data of 
LAGEOS, LAGEOS II and LARES will be finally used by the proponent team in order to try 
to detect the Lense-Thirring effect. This could turn out to be a non-trivial matter because of the 
non-gravitational perturbations. Indeed, if, for instance, a combination 16 

^LARES + ^LAGEOS + /^LAGEOS II ^ 



was adopted instead of that of eq. ( |4. 1| ), the coefficients of the nodes of LAGEOS and LAGEOS II, 
in view of the lower altitude of LARES, would be 

ALAGEOS II ALARES _ ALARES ALAGEOS II 

h, ".2 "4 .2 ".4 IT -391 C 

"1 — ALARES ALAGEOS II _ ALAGEOS IIalaGEoS — 

(4.7) 

ALARES ALAGEOS _ ALAGEOS ALARES 

^ "2 4 2 4 A HHAA 

"2 — ^LAGEOS ^LaGEos H_^LaGEoS ii'a^LAGEOS — 1 

and the combined Lense-Thirring signal would amount to 676.8 mas yr -1 . As a consequence, the 
direct and indirect effects of the non-gravitational 17 perturbations on the nodes of LAGEOS and 
LAGEOS II would be enhanced by such larger coefficients and this may yield a degradation of the 
total obtainable accuracy. 

5. The Sun's Lense-Thirring effect on the inner planets 

The action of the solar gravitomagnetic field on the Mercury's longitude of perihelion 18 UJ = 



£2 + ft) was calculated for the first time by De Sitter in Ref . [ 64 ] who, by assuming a homogenous 



and uniformly rotating Sun, found a secular rate of —0.01 arcseconds per century (arcsec cy in 



the following). This value is also quoted by Soffel in Ref. [14]; Cugusi and Proverbio p2|] yield 



-0.02 arcsec cy 1 for the argument of perihelion ft) of Mercury. 



The recent estimate of the Sun's angular momentum [ |62| , p3Q from helioseismology yields a 
precessional effect one order of magnitude smaller for Mercury; see Table [l4| for the predicted 
Lense-Thirring precessions ET LT of the longitudes of the perihelia of the inner planets; they are 
of the order of 10~ 3 — 10~ 5 arcsec cy -1 . In computing them we have taken into account that the 
inclinations / of the planets which are usually quoted in literature are referred to the mean ecliptic at 
a given epoch 19 [[j~J], while the Sun's equator is tilted by £ = 7.15 deg to the mean ecliptic at J2000 
)8t]. So far, the solar Lense-Thirring effect on the orbits of the inner planets was believed to be too 



I6 The impact of the geopotential is, by construction, unaffected with respect to the combination of eq. (11). 
17 The same may hold also for time-dependent gravitational perturbations affecting the nodes of LAGEOS and LA- 
GEOS II, like the tides. 

18 Since Q. and CO do not lie, in general, in the same plane, tTS is a "dogleg" angle. 
I9 It is J2000 (JD 2451545.0). 
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Table 14: First line: corrections AG7, in 10~ 4 arcsec cy _1 (1 arcsec cy _1 = 10 mas yr = 1.5 x 10~ 15 
s -1 ), to the standard Newton/Einstein secular precessions of the perihelia of the inner planets estimated by 
E.V. Pitjeva [ po| | with the EPM2004 ephemerides. The result for Venus (E.V. Pitjeva, private communica- 
tion, 2008) has been obtained by recently processing radiometric data from Magellan spacecraft with the 
EPM2008 ephemerides [p9[|. The errors in brackets are the formal, statistical ones. Second line: predicted 
Lense-Thirring perihelion precessions GJu, in 10~ 4 arcsec cy~' [100]. Third line: nominal values of the 
perihelion precessions due to the Sun's oblateness for J'f = 2x 10~ 7 [101, the current level of uncer- 
tainty in it is about 10% [102]. Fourth line: nominal values of the perihelion precessions due to the ring of 
the minor asteroids for m^ng = 5 x 10~ 10 M© [ 103 ] ; the uncertainty in it amounts to 8m^ n& = 1 x 10~ 10 M® 
[ 103]. Fifth line: nominal values of the perihelion precessions due to a massive ring modelling the action of 
the Classical Kuiper Belt Objects (CKBOs) for mcKBOs = 0.052 m® jlQ4)]. 
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-4 + 5(1) 
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small to be detected [14]. However, the situation is now favorably changing. Iorio in Ref. [ 106 ] 
preliminarily investigated the possibility of measuring such tiny effects in view of recent important 
developments in the planetary ephemerides generation. 

First attempts to measure the Sun's Lense-Thirring effect have recently been implemented by 



Iorio in Refs. Qlfjg, |107] with the Ephemerides of Planets and the Moon EPM2004 [ |108| ] produced 



by E.V. Pitjeva at the Institute of Applied Astronomy (IAA) of the Russian Academy of Sciences 
(RAS). They are based on a data set of more than 317,000 observations (1913—2003) including 
radiometric measurements of planets and spacecraft, astrometric CCD observations of the outer 
planets and their satellites, and meridian and photographic observations. Such ephemerides were 
constructed by the simultaneous numerical integration of the equations of motion for all planets, 
the Sun, the Moon, 301 largest asteroids, rotations of the Earth and the Moon, including the pertur- 
bations from the solar quadrupolar mass moment jf and asteroid ring that lies in the ecliptic plane 
and consists of the remaining smaller asteroids. In regard to the post-Newtonian dynamics, only 
the gravitoelectric, Schwarzschild-like terms of order ff(c^ 2 ), in the harmonic gauge Q109Q , were 
included; the gravitomagnetic force of the Sun was not modelled. 

The EPM2004 ephemerides were used by Pitjeva in Ref. [^] to phenomenologically estimate 
corrections AtiJ to the known standard Newtonian/Einsteinian secular precessions of the longitudes 
of perihelia of the inner planets as fitted parameters of a particular solution. In Table 14 it is 
possible to find their values obtained by comparing the modelled observations computed using the 
constructed ephemerides with actual observations. In determining such extra-precessions the PPN 
parameters y and j3 and the solar even zonal harmonic coefficient Jf were not fitted; they were 
held fixed to their general relativistic and Newtonian values, i.e. j = f} = I, jf = 2 x 10~ 7 . In 
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the EPM2008 ephemerides Q99Q , used to accurately estimate the perihelion precession of Venus 
as well, were constructed by also modelling the actions of Eris and of the other 20 largest Trans- 
Neptunian Objects (TNOs); the database was enlarged by including, among other things, ranging 
data to Magellan and Cassini spacecraft. 

Although the original purpose 20 of the estimation of the corrections ACT was not the measure- 



ment of the Lense-Thirring effect, the results of Table |14| can be used to take first steps towards an 
observational corroboration of the existence of the solar gravitomagnetic force. 

From Table [R] it turns out that the magnitude of the Lense-Thirring perihelion precessions 
of the inner planets just lies at the edge of the accuracy in determining ACT. All the predicted 
Lense-Thirring precessions are compatible with the estimated corrections ACT. Moreover, if one 
constructs the ratios Jz?ab of the Lense-Thirring precessions 

m A a 3 (1 -e 2 ) 3 / 2 

c/> — LT _ "Bl 1 e BJ sc 

•^ = ^:-4(i-4)3A < 51 > 

for all the pairs of inner planets A and B and compare them with the ratios 

n AB = ^ (5.2) 

of the estimated corrections for the same pairs of planets A and B, it turns out that 

I^ab = =^ab (5-3) 

within the errors 21 . It must be noted that for Venus the formal, statistical error in ACT is smaller 
than the gravitomagnetic effect. Concerning the systematic alias due to the various competing 



dynamical effects listed in Table [14j, the present-day level of mismodelling in them would make 
them not particularly insidious (at « 10% level of accuracy), at least for Mercury and, especially, 
Venus, with the exception of the impact of the Sun's oblateness on Mercury; by the way, it would 
be possible to remove it by suitably designing a linear combination of the perihelia of Mercury 
and Venus, as done for the LAGEOS satellites. In principle, it would be possible to use also 



the nodes [106], if only the corrections AD, to their standard precessions were available. If and 
when other teams of astronomers will independently estimate their own corrections to the standard 
perihelion (and, hopefully, node as well) precessions with different ephemerides, it will be possible 
to fruitfully repeat the present test. 

6. The Lense-Thirring test in the gravitational field of Mars with the Mars Global 
Surveyor spacecraft 



Iorio in Refs. [72, 73] proposed an interpretation of the time series of the RMS orbit overlap 



differences [[71]] of the out-of-plane part v of the orbit of the Martian nearly polar artificial satellite 

2() The goal of Pitjeva in Ref. [ po| was to make a test of the quality of the previously obtained general solution in which 
certain values of jS,/,./^ were obtained. If the construction of the ephemerides was satisfactory, very small 'residual' 
effects due to such parameters should have been found. She writes: "At present, as a test, we can determine [...] the 
corrections to the motions of the planetary perihelia, which allows us to judge whether the values of /3 , y and jf used 
to construct the ephemerides are valid.". The smallness of the extra-perihelion precessions found in her particular test- 
solution is interpreted by Pitjeva as follows: "Table 3 [of Ref. [po|] shows that the parameters y = /3 = l.J^' = 2 x 10~ 7 
used to construct the EPM2004 ephemerides are in excellent agreement with the observations." 

21 In calculating the errors in ITab we linearly propagated the realistic errors in AdJ of Table [l4| not the formal ones. 
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Mars Global Surveyor (MGS) over a time span AP of about 5 years (14 November 1999-14 January 
2005 in J73|]) in terms of the general relativistic gravitomagnetic Lense-Thirring effect. It turned 
out that the average of such a time series over AP, normalized to the predicted Lense-Thirring 
out-of-plane mean shift over the same time span, is jj. = 1.0018 ±0.0053. 



Our interpretation has recently been questioned by Krogh in Ref. [74]. The remarks con 



cerning the analysis presented in Refs. [72, [73|] mainly deal with I) The observable used: Iorio in 
Refs. [[72|, [73|] would have misinterpreted the MGS data II) The confrontation between the predic- 
tion of the gravitomagnetic Lense-Thirring shift and the data over the chosen time span AP: Iorio 



in Refs. [72, 73] would have incorrectly compared the 1.6 m value of the out-of-plane average orbit 



error released by Konopliv et al. in Ref. [|71[] for the entire MGS data set to the Lense-Thirring shift 



calculated for a shorter time interval AP III) The data set used: Iorio in Refs. [72, 73] discarded 



some of the initial months of the MGS data set IV) The bias-neglected by Iorio in Refs. [|72|, [730- 



due to the multipolar expansion of the Newtonian part of the martian gravity field, as pointed out 



in Refs. QllOj , |1 1 1| ] quoted in Ref. [J74J] V) The impact of the atmospheric drag, neglected by Iorio 



in Refs. g 

Below we present our reply which, basically, consists of the following points [j75|]. As fur- 
ther, independent tests, here we present various linear fits to different data sets including, among 
others, the full time series of the entire MGS data (4 February 1999-14 January 2005) as well; 
the predictions of general relativity turn out to be always confirmed. The analytical calculation 
of the competing aliasing effects due to both the gravitational and non-gravitational perturbations, 
which affect the in-plane orbital components of MGS, do not show up in the real data. More- 
over, the non-conservative forces, whose steadily refined modeling mainly improved the in-plane 
orbital components of MGS, not the normal one, exhibit high-frequency, non-cumulative in time 
variations. 

6.1 Our arguments 

1) The entire MGS data set was subdivided by Konopliv et al. in Ref. [jTlJ] in 388 (not 442, 
as claimed by Krogh in Ref. [|74|]) smaller time intervals of data called arcs. For MGS, 
the lengths of the arcs vary from 4 to 6 days, so to cover many orbital revolutions (« 2 h). 
For each arc, the spacecraft position and the velocity, among other things, were estimated 
and used as starting point for a numerical propagation of the satellite's motion by means of 
the dynamical models which, in the case of MGS, did not include the general relativistic 
gravitomagnetic force. Contiguous arcs were overlapped by an amount of just 2 h, i.e. one 
orbital revolution, and the RMS spacecraft position difference among the predicted positions 
propagated from the estimated ones in the previous arc and the estimated positions of the 
subsequent arc was computed. Since the arc overlaps cover just about one orbit, such RMS 
differences, in fact, account for any among measurement errors, random errors, systematic 
bias due to mismodeling/unmodeling dynamical forces yielding secular, i.e. averaged over 
one orbital revolution, effects, whatever their physical origin may be. Indeed, RMS of orbit 
solution overlaps are commonly used in satellite geodesy as useful and significant indicators 



of the overall orbit accuracy [112, 89]. Conversely, they are also used to gain informa- 



tion about systematic errors coming from inaccurate modeling of the forces acting on the 



spacecraft. For details see Refs. [112, 89]. Of course, such a technique is insensitive to 
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short-period effects, i.e. having frequencies higher than the orbital one: only dynamical fea- 
tures of motion with time-scales equal to, or larger than one orbital period can be sensed 
by such orbit overlap differences. Moreover, the average orbit error (AVdiff) of about 1.6 m 
does not refer to this or that particular arc overlap; instead, it comes from the mean of the 
entire set of RMS orbit overlap differences for the chosen time span AP and is well represen- 
tative of those un-modelled/mis-modelled forces yielding effects which do not average out 
over AP, as it is just the case of the Lense-Thirring signal. Time- varying patterns exhibiting 
well-defined periodicities-including also measurement errors like, e.g., those related to the 
Earth-Mars geometry-are, instead, mainly averaged out yielding little or no contribution to 
the average orbit error. Incidentally, from the above discussion about the meaning of the 
average orbit error, it should be apparent that it does not make sense to look for the error of 
the error, as, instead, seemingly required by Krogh in Ref. [[74|] when he blames Iorio [73] 



for not having included the uncertainty in (Av res ). Another criticism in Ref. [ |74| ] is that the 
RMS overlap differences would be unable to specify any orbital precession. 

To reply to all such criticisms we decide to perform another, independent test of our hypoth- 
esis. First, by linearly fitting 22 the full time series of Ref. []7jJ], after having rescaled the data 
points in order to shift the zero point of the time-series to the middle of the data span, we get 
a slope of —0.64 ±0.26 m yr _1 , (with 95% confidence bounds), while the predicted Lense- 
Thirring MGS out-of-plane rate (customarily defined positive along the spacecraft's orbital 
angular momentum) amounts to 0.62 m yr _1 . The obtained minus sign is due to the fact 
that Konopliv et al. in Ref. fl7l| ] defined the normal direction to be positive in the opposite 
direction of the MGS orbital angular momentum (Konopliv 2007, private communication). 
Should such a linear fit be used as indicator of the existence of the Lense-Thirring effect, 
its relativistic prediction would be fully confirmed within the experimental error; instead, 
the hypothesis of a null effect would be rejected at 2.4 sigma level. Then, we also repeat 
our procedure by fitting with a straight line the entire data set without full January 2001, 



mainly affected by likely measurement errors which, according to Krogh y74|], would mimic 
the Lense-Thirring effect, getting —0.61 ±0.26 m yr _1 . The removal of the entire year 2001, 
mainly affected by angular momentum wheel desaturation operations, yields —0.57 ± 0.28 
m yr _1 . Another linear fit to the time series after removing the last month (December 2004- 
January 2005) yields -0.62 ±0.27 m yr" 1 . 

Such results reply to the criticisms II) and III) as well concerning AP, to which a large part 



of Ref. [74] is devoted. 



IV) Krogh [74] quotes Ref. [110] in which analytical calculation about the corrupting impact of 
various physical parameters of Mars through the classical node precessions induced by the 
even zonal harmonic coefficients J( of the multipolar expansion of the Newtonian part of the 
martian gravitational potential are presented. In particular, the authors of Ref. JllOj ] use the 
first five even zonals J2..J\o along with their associated errors from former global solutions 
for the Mars'gravity field, the uncertainty in the Mars' GM and in the MGS semimajor axis 



22 Note that, since the plots in Fig. 3 of Ref. J7l| ] are semi-logarithmic, one should not visually look for a straight 
line in them. 
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and inclination, plug them into analytical formulas for the classical secular node precessions 
and conclude that, since the resulting effect is tens of thousand times larger than the Lense- 
Thirring effect on MGS, this would be fatal for any attempt to detect the gravitomagnetic 
frame-dragging with such a spacecraft. 

The point is that such figures, as others which can be obtained from more accurate calcu- 
lation, must ultimately be compared with the reality of the data, i.e. the RMS orbit overlap 
differences of MGS. 

We, in fact, repeated such calculation by considering also the other even zonals up to /20 
along with the latest errors of the MGS95J global solution and including the uncertainties in 
the Mars' radius as well. By summing, in a root-sum-square fashion, such terms we get a 
mean bias of 78.9 m d _1 in the out-of-plane MGS orbital component: by linearly summing 
them we get an upper bound of 111.6 m d . Such figures clearly show how they are by 
far not representative of the real MGS orbit. Indeed, over a time span of 5 years we would 
have an enormous mean shift as large as 144 km (root-sum-square calculation) or 203 km 
(linear sum). Interestingly, even if the set of the RMS overlap differences of MGS were to be 
considered as representative of a single orbital arc 6 d long only, the conclusion would be the 
same: indeed, in this case, the total cross-track mean shift due to the martian gravitational 
potential would amount to 473.1 m (root-sum-square) or 669.6 m (linear sum). 



In regard to Ref. [Ill], quoted in Ref. [74] as well, let us recall again that the RMS orbit 



overlap differences are just used to account, in general, for all the measuremnt/systematic 



errors giving an indication of the overall orbit accuracy [112, 891]. The important point is that 



they cancel out, by construction, errors, systematic or not, common to consecutive arcs— it 



would just be the case of a bias like that described in Ref. Ql 11| ]— , while effects like the 



Lense-Thirring one, accumulating in time, are, instead, singled out Q89j]. 



V) In regard to the impact of the non-gravitational perturbations, the authors of Ref. Ql 10p yield 
a total un-modelled non-gravitational acceleration of 10~ n m s~ 2 which is the same order 
of magnitude of the Lense-Thirring acceleration induced by Mars on MGS. They neither 
present any detailed calculation of the effect of such an acceleration on the normal portion 
of the MGS orbit nor specify if such a magnitude refers to the out-of-plane component. 
However, some simple considerations can be easily traced: a hypothetic, generic perturbing 
out-of-plane force 6.7 times larger than the Lense-Thirring one and having the same time 
signature, i.e. linear in time, should induce a 10.8 m cross-track shift, on average, over the 
considered time span AP. Again, such a bias is neatly absent from the data. By the way, as 
clearly stated in [[71]], it is the along-track portion of the MGS orbit— left unaffected by the 
Lense-Thirring force— to be mainly perturbed by the non-gravitational forces: indeed, the 
along-track empirical accelerations fitted in Ref. [ f7l| ] amount just to rj 10~ n m s~ 2 , which 



shows that the guess by [110] is somewhat correct, but it refers to the along-track component. 



Time-dependent, periodic signatures would, instead, be averaged out, provided that their 
characteristic time scales are relatively short, as it is just the case. Indeed, the non-conservative 
accelerations, which are especially active in the MGS in-plane orbital components as clearly 



stated in Refs. [|1 13[ [71[], exhibit time- varying patterns over 12 hr [113] which, hypotheti- 
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cally mapped to the out-of -plane direction, are averaged out over multi-year time spans (and, 
incidentally, over 6 d as well). To be more definite, in regard to the issue of the impact of the 



atmospheric drag on the cross-track portion of the orbit of MGS, raised in Ref. [|74]], let us 
note that it requires not only to consider the node £2, as apparently claimed by Krogh J74|], 



but also the inclination i according to Ref. [ 20 ] 



Av 




(A/) 



+ (sin/A£2) ; 



(6.1) 



According to, e.g., Ref. [|30|], the perturbing acceleration Ad rag due to the atmospheric drag 
can be cast into the form ^ 

^drag = -^ZC D — PVV, (6.2) 



where E/M is the spacecraft cross sectional area (perpendicular to the velocity) divided by 
its mass, Cd is the drag coefficient, p is the atmospheric density (assumed to be constant over 
one orbital revolution), v is the satellite velocity in a planetocentric, non-rotating frame of 
reference and Z is a corrective coefficient accounting for the fact that the atmosphere is not 
at rest, but rotates with angular velocity Ya more or less rigidly with the planet; Z 1 for 
polar orbits [fSCQ. While the secular, i.e. averaged over one orbital period T, drag shift on the 



node vanishes, it is not so for the inclination: indeed, it turns out [30] 



(A/) 
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Adrag \ V^A 
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(6.3) 



As a result, the orbital plane tends to approach the planet's equator; the terms in brackets is 
the ratio of the drag force to the Newtonian monopole. As usual in perturbation theory, a is 
meant as evaluated on the unperturbed reference ellipse. Thus, the out-of-plane drag shift is 



from eq. (6.1) 



( A Vdrag> 



rag/ 



V2 



(6.4) 



In the following we will assume that i//a ~ '/'Mars = 7.10x10 5 s . Let us see what happens 
in the (unlikely) worst-case A drag 10~ n m s~ 2 ; it turns out that 



(AVd rai 



g/T 



1 x 10~ 5 m. 



(6.5) 



But A^ag is not constant over time spans days or years long [ 114 ], so that such an effect is 
not a concern here. By the way, even if it was not so, by assuming a 10% mismodeling in 
drag-which is, in fact, modeled in Ref. [f7lj]— eq. (|63|), mapped onto about 5 yr, would give a 
0.7% uncertainty. 



Finally, Krogh [74] remarks that decreasing in the averages of the RMS orbit overlaps oc- 



curred in view of constantly improved modeling [ |113| , |1 15J , but he does not recognize that 
the better modeling of the non-gravitational forces acting on MGS introduced in Ref. 



with respect to previous works [113, 115] in which the Lense-Thirring effect was not mod- 
elled as well, only affected in a relevant way just the along-track RMS overlap differences 



(a factor 10 better than in Refs. [113, 115]), not the normal ones (just a factor 2 better than 
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in Refs. [113, 1 15]). Moreover, if the relativistic signature was removed or not present at all 
so that the determined out-of-plane RMS overlap differences were only (or mainly) due to 
other causes like mismodeling or unmodeling in the non-gravitational forces, it is difficult 
to understand why the along-track RMS overlap differences (middle panel of Figure 3 of 
Ref. [f7l|]) have almost the same magnitude, since the along-track component of the MGS or- 
bit is much more affected by the non-gravitational accelerations (e.g. the atmospheric drag) 
than the out-of-plane one. 



7. Conclusions 



• The present and future tests of the Lense-Thirring effect in the gravitational field of the Earth 
with some existing and future SLR artificial satellites are challenged by several competing 
classical forces whose impact is difficult to be assessed in a realistic and conservative way. 

- Because of the lingering statistically significant discrepancies among the estimated val- 
ues of the low-degree even zonal harmonics of the geopotential in several GRACE- 
based global solutions produced by different institutions worldwide, the present-day 
uncertainty in the Lense-Thirring test performed by linearly combining the nodes of 
the LAGEOS and LAGEOS II terrestrial satellites is larger than the claimed value of 
5 — 10% by a factor as large as up to 2-3 times. 

- A different methodology to extract the Lense-Thirring signal from the data of the 
LAGEOS -type satellites should be used to make a complementary test which would 
enforce (or disprove) those performed so far. That is, the gravitomagnetic force should 
be explicitly included in the dynamical force models of LAGEOS and LAGEOS II, and 
a dedicated parameter accounting for it should be estimated in the least-square fitting of 
the observations of the LAGEOS satellites; moreover, the changes, if any, in the values 
of the other estimated parameters with respect to the case in which the gravitomagnetic 
force was not modelled may be inspected as well. Alternatively, corrections AQ. to 
the standard Newtonian node precessions of the LAGEOS satellites may be estimated 
in a purely phenomenological way without modelling the Lense-Thirring effect; their 
values should, then, be compared with the predicted gravitomagnetic node precessions. 
An analogous procedure was followed with the periastron rates in the binary pulsar 
systems and with the perihelion precessions of the inner planets of the Solar System, 
although for different scopes. Finally, another feasible approach, at least in principle, 
consists of producing a new global gravity solution by re-processing the GRACE ob- 
servations with dynamical force models including the gravitomagnetic force as well, 
and estimating it along with all the even zonals of the geopotential and inspecting the 
resulting covariance matrix. 

- The LAGEOS-type LARES satellite, approved by ASI, should be launched at the end of 
2009 or at the beginning of 2010 with a VEGA rocket: its goal is a measurement of the 
Lense-Thirring effect with an accuracy of the order of 1%. Unfortunately, its orbital 
configuration will be different from the originally proposed one, with an altitude of 
about 1400 km; this may pose serious problems in terms of the impact of the even zonal 
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harmonics of the geopotential because a much larger number of them should come into 
play by degrading the total accuracy obtainable. Extensive calculations performed with 
standard geodetic techniques and several global gravity model solutions point towards 
a systematic uncertainty which may be orders of magnitude larger than the claimed 1%. 
Also certain subtle issues related to the indirect impact of the atmospheric drag on the 
node precessions may contribute to further increasing the total error. 

• Recent improvements in the field of the planetary ephemerides make the possibility of mea- 
suring the Lense-Thirring effect in the gravitational field of the Sun with the inner planets 
much more promising than in the past. At present, the magnitude of the Lense-Thirring 
perihelion precessions, not included in the models of the dynamical theories with which the 
planetary observations are reduced, lies just at the edge of the accuracy in determining the 
secular perihelion precessions of the rocky planets. It turns out that the predicted gravito- 
magnetic rates are in agreement with the recently estimated extra-precessions of the perihelia 
for all the inner planets; the same also holds for the ratios of the perihelia for all the pairs 
of inner planets. Concerning the systematic errors which may be caused by the imperfect 
knowledge of several competing forces of classical origin, the most insidious one is the 
Sun's oblateness. However, its impact can be removed by linearly combining the perihelia 
of the planets as done for the LAGEOS-LAGEOS II test. The best suited candidates are 
Mercury and Venus whose orbits should be determined with an higher level of accuracy in 
the near future when ranging data to present and future planned spacecraft orbiting them will 
be collected and processed. It would also be important if other teams of astronomers inde- 
pendently would estimate their own corrections ACT to the usual perihelion (and, hopefully, 
node as well) precessions of the inner planets with different ephemerides to repeat the tests 
discussed here. 

• Another Solar System scenario which could also be used for testing the Lense-Thirring effect 
is the gravitational field of Mars. In this framework, we recently proposed an interpretation 
of some orbital data of the nearly polar Mars Global Surveyor probe in terms of the gravito- 
magnetic force of Mars by performing two independent preliminary tests. The average (Av) 
of the time series of the out-of-plane RMS orbit overlap differences of MGS, normalized 
to the average Lense-Thirring shift (Avlt) over the same time span, is }i = 1.002 ±0.005. 
Linear fits of the time series for different time spans, i.e. including different sets of data 
points, yield results in agreement with the LT interpretation. For example, the fit of the en- 
tire data set yields a slope of 1.03 ±0.42, in LT-normalized units; removing the data for the 
month January 2001, likely affected by major measurement errors, yields a normalized slope 
0.98 ±0.42. 
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